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Multivariable generalizations of the continuous Hahn and Wilson polynomials are introduced as 
eigenfunctions of rational Ruijsenaars type difference systems with an external field. 

PACS numbers: 02.30.Gp, 03.65.-w, 02.90.+p 

In a by now famous pioneering paper, Calogero studied (essentially) the spectrum and eigenfunctions of a quantum 
system of N one-dimensional particles placed in a harmonic well and interacting by means of an inverse square 
pair potential Q. The spectrum of the system, which is discrete due to the harmonic confinement, turns out to be 
remarkably simple: it coincides with that of noninteracting particles in a harmonic well up to an overall shift of the 
energy. The structure of the corresponding eigenfunctions is also quite simple: they are the product of a factorized 
(Jastrow type) ground state wave function and certain symmetric polynomials. Recently, a rather explicit construction 
^ , of these polynomials was given in terms of raising and lowering operators ||. 

Some time after its introduction, Olshanetsky and Perelomov realized that the Calogero system can be naturally 
^\ , generalized within a Lie-theoretic setting, such that for each (normalized) root system there exists an associated 
t-H ■ Calogero type quantum model ||. From this viewpoint, the original model corresponds to the root system Ajv_i. 
If one restricts attention to classical (i.e., non exceptional) normalized root systems, then the Hamiltonians of the 
corresponding Calogero models split up in two types: 
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(dj — d/dxj). For the very special case N = 1, (the polynomial parts of) the eigenfunctions of these Hamiltonians 
amount to Hermite polynomials (type A) and Laguerre polynomials (type B). If the external harmonic field is switched 
off, i.e. for ujq = 0, then the spectrum of the Hamiltonian becomes continuous (assuming g,go > 0), and the 
eigenfunctions are no longer polynomials but give rise to multivariable families of Bessel type functions B . 

Several years ago, Ruijsenaars introduced a relativistic generalization of the type A model without external field 
The Hamiltonian of this relativistic system is given by an analytic difference operator that goes over in Hcai.A 
(|l|) (with loq — 0) after sending the step size to zero. (This transition may be interpreted as the nonrelativistic limit.) 
More recently we found similar difference versions for the type A model with cjq ^ and for the type B model [[|. 
In this letter we study the eigenfunctions of these two difference models. For the type A difference version this will 
lead to a multivariable generalization of the continuous Hahn polynomials [QH , whereas for type B we will find a 
multivariable version of the Wilson polynomials |^|,[l0|. By sending the step size to zero, we recover the spectrum and 
eigenfunctions of the Calogero Hamiltonians Hcai.A (P and Hc a i,B 
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Type A: continuous Hahn case. The Hamiltonian of the (rational) difference Ruijsenaars system with external 
field is given by the second order difference operator |j] 

1< 3 <N k^j kjij 

+w 1 !\x J ) H v l ' 2 {x k - Xj ) e id i J] v^ixj - x^fix,) 

-w+(xj) Y[v(xj -x k ) - w-(xj) Y[v(x k - xj)) (3) 

kjij k^j 

where (e ±ldj ^>)(xi, . . . , ,xn) = *&(xi, . . . ,Xj-i,Xj ± i,xj+i, . . . ,xn) and 

v(z) = 1 + g/(iz), w + {z) = {a + + iz)(b + + iz), W-{z) = (a_ — iz)(b- — iz). (4) 

We will assume 

g>0, a_=a + , 6_=6 + , Re(a + , a_, b + , > 0, (5) 

which ensures, in particular, that the Hamiltonian is formally self-adjoint. In order to solve the eigenvalue problem 
for H A (§)-(§) in the Hilbert space of square integrable permutation invariant functions, we introduce the weight 
function 

A A = [] T (9 + ^-x k )) r(fl+ + .^ )r(5+ + ix . )r(fl _ _ lXj)T{b _ _ ix .y (6 ) 

l<jjtk<N [l[Xj Xk)) l<j<N 

Condition (|^) implies that is positive. Since the transformed operator 

Ha = A- 1/2 H A A\ /2 

= ( w +( x i) II v ( x 3 ~ x k) ( e ~' ld3 - !) + 

l<j<N k^j 

w^{x J )\{v{x k -x ] ){e^ -1)) (7) 
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clearly annihilates constant functions, it follows that is an eigenfunction of H A with eigenvalue zero. This 

1 /2 

eigenfunction corresponds to the ground state. The excited states are a product of A^ and symmetric polynomials 
associated with the weight function A^ (^|). 
Specifically, one has 

H A ^n = E A {n)^ft, n<EZ N , nx>n 2 >--->n N >0, (8) 

with eigenvalues reading 

E A (n)= Y nj(nj + a+ + a- + b + + b--l + 2(N-j)g), (9) 

l<j<N 



and eigenfunctions of the form 



= A A /2 PnyA {x), (10) 



where Pu,a(x) denotes the symmetric polynomial determined by the conditions 
A.l Pa, A (x) = mfi(x) + ^ c «,«' m n>(x), Cfifi' £ C; 



A. 2 / pn yA (x) irifii(x) A A dx\ ■ ■ ■ dx^ = if ft' < fi. 
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Here the functions m,f L (x) denote the basis of monomial symmetric functions 

nef, m > n 2 > ■ ■ ■ > n N > 0, 



TV ' 



and the partial order of the basis elements is defined by 



n < n iff n'j < nj for k = 1, . . . , N 

l<j<k l<j<k 



(11) 



(12) 



(ft' < ft if ft' < ft and ft' ^ ft). 

The proof of the above statement boils down to showing that the polynomials pfi^(x) are eigenfunctions of the 
transformed operator Ha @)- This follows from the fact that Ha is both triangular with respect to the monomial 
basis: 



(13) 



and symmetric with respect to the L 2 inner product with weight function (^). The eigenvalues EA(n) @ 
are obtained by computing the diagonal matrix elements \Ha]u,r, i-e., the leading coefficients in Expansion (113) of 
(7iA TO ra)(^ f ) m monomial symmetric functions mff(x). 

The polynomial pn,A(x) by definition amounts to mft(x) minus its orthogonal projection in L 2 (R N , A^cfei • • • dxw) 
onto span{m,-j'}fj' <? -j. For N = 1 these polynomials reduce to (monic) continuous Hahn polynomials 

Type B: Wilson case. Multivariate Wilson polynomials are obtained much in the same manner as their contin- 
uous Hahn counterparts, except that in addition to being permutation symmetric now everything also becomes even 
in Xj, j = 1, . . . , N. The Hamiltonian of the type B version of the difference Ruijsenaars system reads || 
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and pfi^six 2 ) is the even symmetric polynomial determined by the conditions 
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B.l p,xb(x 2 ) = m rl (x 2 ) + ^2 c a,n> m R ,(x 2 ), c n . n > G C; 

n' <n 

B.2 / pff b(x 2 ) mft'ix 2 ) Ab dxi ■ ■ ■ dxiy = if n' < n. 
Here the functions nif L (x 2 ) stand for the basis of even symmetric monomials 

m R (x 2 ) = xf 1 ' 1 ■ ■ ■ x 2 "' N , ne1 N , n 1 >n 2 >--->n N >0, (21) 

«'GSjv(n) 



and the partial order of the basis elements is the same as before (see (12)). 

The polynomials ph,b(x 2 ) are of course again eigenfunctions of the transformed operator 

Kb = A B 1/2 H B A B /2 



l<j<N k^j 

w (~Xj) n v(-Xj + x k ) v(- Xj - x k ) {e id * - 1)) . (22) 



For N = 1 they reduce to (monic) Wilson polynomials lOUlCf] . 

Transition to the Calogero system. If we substitute Xj — > (3~~ 1 Xj (so 9j — > and a+,a_ — > (/3 2 a;o) — 1 , 

6_l_, fo_ — > (/3 2 Wn) _1 in our type A difference Hamiltonian and multiply by f3 u>oU)' , then we arrive at a Hamiltonian of 
the form Ha (||) with exp(±idj) replaced by exp(±z/39j) and 

v(z) = l + 0g/(iz), w±(z) = /3- 2 (l±if3uJo z)(l ±1(3^ z). (23) 

By sending the step size (3 to zero the difference operator goes over in a differential operator of the form Hoai^A — £q,a-, 
where Hc a l,A is given by ([!]) with coo replaced by coo + u^, and £q,a is a constant with value 

e , A = (uJo + u/ )N{l + (N- l)g). (24) 

In this limit the weight function (^) goes (after dividing by a divergent numerical factor) over in 

^Cal,A= 1] l%-^l 9 II e-^+^ x \ (25) 

Thus, we recover the spectrum and eigenfunctions of the type A Calogero system: 

Hc a i,A i &n = Ecai,A( n )*n, ft G Z N , n x > n 2 > ■ ■ ■ > n N > 0, (26) 

with 

E C ai,A{n) = sq,a + 2(w + lo' q ) ^ n h ^n(x) = A]ll l A pn,Cai,A(x), (27) 

l<j<N 

where Pn Cai,A(x) denotes the symmetric polynomial determined by Conditions A.l, A. 2 with A^ ^ replaced by 

A C al,A 

For type B, the transition to the Calogero system is quite similar. The substitution Xj — ► fi~ x Xj (dj — > /3d j), 
a — > go; & - * .9o + 1/2, c — > (/3 2 wo) _1 : d — v (/3 2 o;q) — 1 , and — > AP 2 ujooJqHb leads to a Hamiltonian of the form 
( |l4| ) with exp(±i<9j) replaced by exp(±i/3<9j) and 

= l + pg/(iz), 

w{z) = /T 2 (l + ^)(1 + ( .^ g ° /2) )(l + ifao z)(l + */J</ z). (28) 

For /3 — > one now obtains a differential operator of the form Hc a i,B — £q.b, where Hcai.B is given by (0) with go 
replaced by g + g' and lo q by oj + w'q, and the constant e 0i s reads 
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£o,B = (w + u/ )N(l + 2(N - l)g + 2(g + g' Q )). 



(29) 



In the limit (and after dividing by a divergent numerical factor) the type B weight function goes over in 

A C aLB = J] IXj+X^lXj-Xkl" J] \ Xj | 2 ^°+So) e^^o)^ (30) 
i<j^k<N l<j<N 

and the eigenfunctions become: 

H Ca l,B^H = E Ca l,B(n)^H, n e Z N , n x > n 2 > ■ ■ ■ > n N > 0, (31) 

with 

Ecai,B(n) = e ,s + 4(w + uj' ) nj, = B Pn,Cai,B(x 2 ), (32) 

l<j<N 

where Pft,Cai,B(x 2 ) denotes the even symmetric polynomial determined by the Conditions B.l, B.2 with As j2tj| ) 
replaced by A Ca i,B ©■ 

Let us conclude by remarking that both our multivariable continuous Hahn and Wilson polynomials are limiting 
cases of a multivariable version of the Askey- Wilson polynomials (T^] introduced by Koornwinder Q] as a gener- 
alization of Macdonald's polynomials associated with the root system BCn (f^Q ■ Koornwinder 's multivariable 
Askey- Wilson polynomials are joint eigenfunctions of an algebra of commuting difference operators with trigono- 
metric coefficients ]l5| . This algebra constitutes a complete set of quantum integrals for a difference version of the 
trigonometric BCV-type Calogero-Sutherland system. Similar algebras consisting of commuting difference operators 
that are simultaneously diagonalized by our multivariable continuous Hahn and Wilson polynomials can be obtained 
as rational degenerations. The difference Hamiltonians considered in this letter are the simplest (i.e., lowest order) 
nontrivial operators in these algebras. 
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